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Abstract 
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BB pairs. The branching-fraction upper limits are set at }3{B~^ — > e^u^) < 7.9 x 10^^ and 
B{B+ fi+iy^) < 6.2 X 10~^ at the 90% confidence level. 
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1 INTRODUCTION 



Leptonic decays of B mesons are of interest both for their use in characterizing Standard Model 
(SM) processes and as probes of possible physics beyond the SM. They are, however, experimentally 
challenging to observe due to both their small branching fractions and the difficulty of measuring 
decay modes with neutrinos in the final state. In the SM, these leptonic decays may proceed via 
tree-level processes mediated by a virtual boson. The SM branching fraction for this type of 
decay [1] is given by 



where Gf is the Fermi coupling constant, m; and niB are the lepton and S-meson masses, and 
tb is the B lifetime. Measurements of the leptonic branching fractions can therefore constrain the 
product of \Vui,\ and /b, where \Vub\ is the Cabibbo-Kobayashi-Maskawa matrix element parame- 
terizing the b ^ u coupling and Jb is a decay constant describing the wavefunction overlap of the 
two quarks of the B meson. has been obtained from analyses of semileptonic B decay modes 
[2]. A measurement of the branching fractions in Eq. 1 would allow us to experimentally constrain 
/b, which can be calculated in lattice QCD with uncertainties of ~15% [3]. Measuring a prod- 
uct \Vub\ X /b, that is inconsistent with measurements using other methods, could be interpreted 
as evidence for new physics. An enhancement to the branching fraction could be caused by the 
existence of additional contributions, for example as predicted by supcrsymmetry [4]. 

As shown in Eq. 1, the decay rates arc hclicity suppressed by a factor of mf/m^, yielding 
SM branching fraction predictions of C(10^^^) and 0(10^'') for the e and modes respectively. 
The most stringent upper limits on the leptonic decay modes B~^ — > e"*" Ve and B~^ — > /j,'^ are 
1.5 X 10~^ [5] and 6.6 x 10~^ [6], respectively, at the 90% confidence level. A preliminary results of 
B{B+ e+Ue) < 5.4 x 10"^ and B{B+ fi+u,^) < 2.0 x lO"'^ at the 90% confidence level [7] are 
also available from the Belle Collaboration. 

We search for the decays B'^ e^v^ and B^ IJ-'^^n using a technique in which the accom- 
panying B~ is reconstructed exclusively in one of several hadronic decay modes. This technique 
has previously been applied to other rare-decay searches [8, 9], but it has not been used previously 
in a search for B^ i'^i^e- Although the hadronic-reconstruction procedure has a relatively low 
efficiency, this method has the advantages of highly suppressed backgrounds and knowledge of the 
signal-lepton energy. Because it is very unlikely to observe any background events, an observation 
of events in the signal region would be highly suggestive of a signal. 

The leptonic decay B^ — > t^Ut is expected to have a branching fraction of O(10~^) [10]. 
Although this decay is less suppressed than the e and /x modes, analysis of it is more difficult due 
to additional neutrinos in the final state. The preliminary result from Belle Collaboration (using 
a similar tagged- i? method) reports first evidence for B~^ — > t^z^t, with a branching fraction of 
1.061 

ais (stat) loj6 (syst) xlO [11]. As the e and modes are experimentally cleaner, they may 
ultimately yield more accurate branching-fraction results than S+ — > t'^Vt even with the higher 
suppression. 



The on-resonance data sample used in this analysis corresponds to an integrated luminosity of 
208.7 fb~^, accumulated at the T{AS) resonance. The events were recorded by the BABAR detector 




(1) 



2 THE BABAR DETECTOR AND DATASET 
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at the PEP-II asymmetric e^e~ storage ring. In addition to the on-resonance sample, an ofF- 
resonance sample of 21.5 fb"''' was recorded approximately 40MeV below the T(4S) resonance, 
which is used for studies of backgrounds not originating from the decay of T{AS). 

The BABAR detector is described in greater detail in the literature [12]. Charged-particle track- 
ing and dE/dx measurements for particle identification are provided by both a five-layer double- 
sided silicon vertex tracker and a 40-layer drift chamber contained within the magnetic field of a 
1.5 T superconducting solenoid. A ring- imaging Cherenkov detector (DIRC) provides efficient parti- 
cle identification. The energies of neutral particles are measured by an electromagnetic calorimeter 
consisting of 6580 CsI(Tl) crystals. Muon identification is provided by resistive plate chambers. 
Signal efficiencies and background rates are estimated using a Monte Carlo (MC) simulation of the 
BABAR detector based on GEANT4 [13]. 

3 ANALYSIS METHOD 

The decay B"*" — ^ t-'^i'i produces a mono-energetic charged e or /x in the B"*" rest frame, accompanied 
by the missing-energy signature of a neutrino. To gain sensitivity to this signature, a B meson from 

the decay T(4S') BB, Btag, is reconstructed in hadronic decay modes B"*" D^*^^X~^, B~ — > 
BO B)(*)-X+ and B^^ D^*^+X- [14]. While both charged and neutral Btag candidates 
are considered for reconstruction, the events where the best Btag candidate (as discussed below) is 
neutral are vetoed. B° candidates are reconstructed in the modes B° iC~7r+, B° K^tt+tt^, 
Z)0 ^ K- '7r"'"7r"'"7r and ii'g7r"'"7r ; D*^ candidates are formed from a and either a vr'^ 

or a 7. Similar reconstruction is done for the D^*^^ candidate. Pions and kaons remaining in the 
event are combined in subsets with the D candidate to form a Btag candidate, considering the 
quantity |AB| = \Eb — Bbeamj, where Eb is the energy of the reconstructed B and Bbeam is the 
beam energy, both in the ccntcr-of-mass (CM) frame. Btag candidates must satisfy the requirement 
I AB| < 0.2 GcV, and the Btag with the lowest value of | AB| is chosen as the candidate. In the case 
of multiple Btag candidates, the one with a highest a priori purity is chosen. The energy-substituted 
mass 

mES = ^JeI.^-p'b, (2) 

where pB is the momentum of the B candidate in the CM frame, is required to be 5.270 GeV/c^ < 
n^ES < 5.288 GeV/c^. The tues distribution is shown in Fig. 1. The resolution in tubs is dominated 
by the spread in Bbeam, ^-Ebeam ~ 2.59 MeV. 

Event shape requirements are used to suppress the combinatorial backgrounds not originating 
from the decay of T{4S). These requirements are: R2 < 0.5, where R2 is the ratio of zeroth and 
the second Fox- Wolfram moment [15]; ] cos 0^1 < 0.9, where 9t represents the angle between the 
thrust axis defined by the Btag candidate and by the combination of all other tracks and clusters 
in the event. 

The tracks and clusters not used in the Btag reconstruction are assumed to originate from the 
decay of the signal B, Bgignai- Since the CM energy is precisely known, reconstruction of the Btag 
fully determines the energy and momentum, and thus the rest frame, of the Bgig^al candidate. 
The signal lepton energy is constrained to be at the kinematic endpoint in the Bgignai rest frame, 
producing a distinct experimental signature. Fig. 2 shows a comparison of the lepton-momentum 
distributions in the Bgignai rest frame and the CM frame. The Bgignai rest frame is calculated from 
the recorded beam energies, and a momentum opposite in direction to the reconstructed Btag- We 
require the highest momentum track of correct charge in the Bgignai frame to have a momentum 
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Figure 1: mEs distribution of -Btag candidates. Signal MC events are shown on top and background 
MC events on the bottom. The events are required to have one reconstructed charged B. The signal 
and background samples are scaled to data luminosity, assuming both the signal mode branching 
fractions to be 2.0 x 10~^ (i.e., the current upper limit of the mode — > n'^v^). 



p* between 2.54 and 2.72 GeV/c^. This highest momentum track is also required to satisfy (e or 
/i) lepton identification requirements, and to survive a veto on kaon particle identification based 
primarity on DIRC. 

The missing momentum in the i^signai frame is measured to be 

P'^miss = P*T(iS) - P*Bt.s ~ P*i ' (3) 

i 

where the sum is of all charged and neutral particles not associated with the reconstruction of 
the i?tag) and all quantities are computed in the -Bgignai rest frame. The direction of the missing 
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Figure 2: The lepton-candidate momentum distributions in the -Bsignai frame (top) and in the CM 
frame (bottom) before applying the signal-selection criteria. The resolution gain provided by the 
Btag reconstruction allows a tighter selection around the signal peak. The signal MC is scaled 
arbitrarily. 



momentum is described by 

CO«^Pmi.. = f^, (4) 
IPmiss I 

where Pz^niss component of the momentum relative to the axis defined by the e~ (high 

energy) beam direction, and the momenta are computed in the CM frame. To ensure that the 
missing momentum is carried by the neutrino rather than a particle passing outside of the detector 
acceptance, we require —0.76 < cos6'p^j^^ < 0.92. 

We also require that the missing momentum be consistent with a neutrino recoiling against the 
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signal candidate lepton by requiring AP^j^^ < 0.7GeV/c, where 

^-^miss IPmiss ~^ P I • 



(5) 



This requirement is useful in suppressing semileptonic B decays, where a high momentum lepton 
may be present. 

The amount of extra energy in the event is described by a quantity ^^extraj defined by 



where Ei is the CM frame energy of a given particle, and the sum runs over all tracks and clusters 
not associated with the Stag or the lepton candidate. Wc select events with i^cxtra < 1.2 GeV. 
The resulting selection efficiencies after applying the above selection criteria are (0.122 it 0.012) % 
and (0.145 ± 0.013) % for e and ^ modes, respectively, where the uncertainties include only MC 
statistics. 

Backgrounds potentially arise cither from sources with a peak in Stag n^'ESi or which are 
combinatorial in -Stag n^'ES- One source of the nonpeaking background is from the misreconstruction 
of the decay products of r+r~ or qq {q = u,d,s,c). These backgrounds are strongly suppressed 
by event shape requirements. Nonpeaking background can also arise from misreconstructed Btag 
candidates in BE events. 

No nonpeaking background events in MC (t~^t~ and qq samples are approximately equal to 
the data luminosity, and BB sample is ~ 5 times the data sample) pass the selection criteria. The 
absence of nonpeaking background is verified by studying the p* and niES sidebands (2.0 GeV < 
p* < 2.5GcV and 5.22GcV/c^ < m,ES < 5.26GeV/c^). The ratios of the number of events in these 
sidebands are used to extrapolate the background level in the signal region. For B^B~ events, 
the nonpeaking contribution is estimated by fitting the ruES distribution with a combination of a 
Gaussian and an ARGUS function [16]. The background estimate from the sidebands is consistent 
with the prediction from the MC of no nonpeaking background. As can be seen in Fig. 1, the 
niES sideband region is well modeled by background MC up to the uncertainties arising from the 
relative cross sections of various decay modes. For nonpeaking events, the quantities p* and niES 
are assumed to be uncorrelated in both the signal and sideband regions. 

In peaking background events, Stag is correctly reconstructed, but the B^\gaai decays via a 
non-signal mode. The most likely backgrounds are h — > uiv^ events and two-body decays with a 
misidentified tt, where the signal lepton may be near the kinematic endpoint. The total B^B~ 
background (nonpeaking and peaking together) is estimated by counting the number of events in 
the MC sample that pass the selection criteria. Requirements on the signal-lepton momentum and 
the missing momentum are especially important in discerning the signal events. After rescaling 
as discussed below, the B'^B~ contribution to the background is estimated to be 0.1151q o75 and 
0.2291q;]^42 events per 209 fb~^, for e and /x modes respectively. Since the contributions from the 
other background types are found to be negligible, the above quantities are taken as the total 
background estimates. 



The study of the lepton-momentum sidebands reveals a slight discrepancy between the niES peak- 
ing B^B~ yields in data and in MC. To estimate the Stag yield in the samples, the ruES distribution 




(6) 



4 SYSTEMATIC STUDIES 
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for the on-resonance data sample is studied after subtracting the normahzed off-resonance rriES dis- 
tribution from it. To compensate for the CM energy difference, the off-resonance niES distribution 
is shifted by ~ 20MeV, matching the kinematic endpoints of the two distributions. The precise 
vahies of the required energy shift are determined by comparing the on-resonance and off-resonance 
CM energy distributions. By fitting a sum of a Gaussian and an ARGUS function to the back- 
ground subtracted tues distribution, we determine the parameters for the shape of the nonpeaking 
contribution. Another set of such parameters is estimated by a similar fit using a simulated BB 
sample. These parameters are used to fit the signal region in the data and MC BB samples. The 
yields are calculated by integrating the Gaussian component of the resulting fits, resulting in a 
range of values for the yield (see Table 1). The scaling factor, which is the ratio of the data and the 
MC yields, is calculated to be 0.908 it 0.013{stat) it 0.015{syst), where the systematic error arises 
from the uncertainties involved in the fitting procedure. It is notable that while the yield itself 
has a large uncertainty of ^ 8.8%, the ratio between the data and MC yields are less dependent 
on the exact shape of the nonpeaking contribution. Applying the correction factor to quantities 
determined from the BB MC sample introduces a systematic error. However, the signal efficiency 
and the background estimates determined from MC samples are statistically limited. 

Table 1: The peaking BB yields resulting from using different ARGUS function shapes for describ- 
ing the combinatorial contribution. The yields are expressed in B~^B~ pairs per fb~^. 



ARGUS shape source 


MC yield 


Data yield 


correction 


From Data 


2253.38 


2018.25 


0.896 


Prom MC sample 


2581.89 


2380.38 


0.922 


Average 






0.908±0.013 



Tracking efficiency uncertainties are considered, introducing an additional 0.8%/track system- 
atic error on the signal efficiency and the background estimate. The lepton-identification efficiencies 
for data and MC are compared and correction factors of 0.96ib0.01 and 0.962ib0.015 are determined 
for e and n respectively. This correction factor affects all quantities determined from MC samples, 
introducing a systematic uncertainty. The misidentification rate of leptons as pions is also studied. 
The difference between the misidentification rates between data and MC introduces a correction 
factor to background estimates. In the case of muons, the pion misidentification rate is estimated 
to be (5.0 lb 0.5) % for MC and (5 it 1) % for data. For the electron identification, the discrepancy 
is larger with misidentification rates of (0.01 ±0.01) % in MC and (0.05 ±0.02) % in data. This dif- 
ference requires the number of background events, estimated from MC, to be scaled appropriately. 
However, as all MC events that pass the selection criteria are verified to have a correct lepton 
identification, we do not include a misidentification-rate correction to the background estimate. 
Additional systematic errors, associated with the quantities -Egxtra and APj^^^g, are estimated by 
varying the selection criteria by 100 MeV. The variations in the signal efficiencies are found to be 
very small, requiring no further additions to the systematics. 

The efficiencies and the background estimates, after applying all corrections, are listed in Ta- 
ble 2. The etag and Cgig correspond to the efficiencies of the i?tag reconstruction in the signal region 
and the signal-lepton selection respectively. The etot represents the product of these two quantities. 
The quantities iVbg; -^SM and iVobs represent, respectively, the number of estimated background 
events, the expected number of signal events according to the SM branching fractions and the 



13 



number of signal candidates observed in the data. 

Table 2: Efficiencies and background estimates, used for calculating the branching fractions, after 
applying all corrections. The first error is statistical, the second is systematic. 



Quantity 


5 ^ /i+ l/^ 


B ^ e+ Ue 




l.OSnb 




208.7/6-1 




229.953x10'' 


etag (%) 


0.239 ±0.013 ±0.004 


0.247 ±0.013 ±0.004 


esig (%) 


60.5 ±4.0 ± 1.0 


49.4 ±3.8 ±0.8 


etot (%) 


0.145 ±0.013 ±0.003 


0.122 ±0.012 ±0.003 


iVbg 


0.229t|]:}^^ ± 0.007 
~ 0.03 


0.115+|]:J^^ ±0.004 
~ 3 X lO"'^ 


iVobs 









5 PHYSICS RESULTS 

Fig. 3 shows the momentum distribution of the signal lepton in the Bsignai frame for data and 
background MC. The absence of events in the signal region is consistent with SM expectations. 




Lepton momentum (B frame) (GeV/c) 

signal 

Figure 3: Signal-lepton momentum in -Bsignai frame after all signal selection criteria. No events in 
data are present in the signal region. 



The branching fraction B is given by 



B{B^ - t.,) = , (7) 
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where A^obs is the number of events that pass the selection, A'^bg is the estimated background count, 
N^+B- is the number of B^B^ pairs in the data sample, and etot is the total efficiency of the 
selection given by signal MC. The upper limit on B at the 90% confidence level is determined 
using a frcquentist procedure that incorporates systematic uncertainties in the signal efficiency and 
expected number of background events [17]. This procedure yields B{B~^ — > e+fg) < 7.9 x 10~^ 
and B{B+ /x+f^) < 6.2 x 10"^ at the 90% confidence level. 

6 SUMMARY 

We have set the branching fraction upper limits for rare leptonic decays B^ —>■ and B~^ 

/x+ to be 7.9 x 10~^ and 6.2 x 10~^ respectively, at the 90% confidence level. Both results are 
consistent with the Standard Model and represent improvements to the most stringent published 
upper limits to date. 
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A Appendix 

Appended are additional plots of interest. 
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Figure 4: The distribution of the lepton candidate momentum in the -Bgignai frame. The events are 
required to pass all reconstruction cuts. The signal —>■ fi'^ Ufj_ and the backgrounds are scaled 
to onpeak data luminosity, assuming SM predictions for branching fractions. The signal — > e"*" 
Ue is scaled with respect to — >■ by the relative luminosity of the sample size. 
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Figure 5: The distribution of AP^^^^, for signal and background samples. The events are required 
to pass all reconstruction and signal cuts (with a relaxed requirement of p* > 2.0 GeV to increase 
background statistics). The signal B'^ i/^ and the backgrounds are scaled to onpeak data, 

assuming SM predictions for branching fractions. The signal — > e"*" I'e is scaled with respect to 
B~^ —I- /i+ i/^ by the relative luminosity of the sample size. 
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Eextra (CM frame) (GeV) 

Figure 6: -Eextra, the total energy not accounted for by the signal lepton candidate, distribution. 
The events are required to pass all reconstruction and signal cuts (with a relaxed requirement of 
p* > 2.0 GeV to increase background statistics). The signal — > //+ and the backgrounds are 
scaled to onpeak data luminosity, assuming SM predictions for branching fractions. The signal 
e"*" fe is scaled with respect to B~^ /x"*" by the relative luminosity of the sample size. 



19 



